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Abstract 

We study the dimuon yield from open charm decays in the intermediate mass 
region(IMR) in heavy-ion collisions at the CERN-SPS energy. We find that 
final state rescatterings which broaden the rriT spectrum of charmed mesons 
can enrich the part of the phase space covered by the NA50 experiment, thus 
leading to an apparent enhancement of IMR dileptons within the acceptance 
of the experiment. Such an enhancement increases with the dimuon invariant 
mass and the single muon energy cut-off. 


Dileptons are important signals in high-energy heavy-ion collisions. Depending on their 
invariant masses they can probe the collision processes at different stages. Low invariant- 
mass dileptons are predominantly produced at the late stage of the collisions. So they 
provide information about the medium modihcation of hadron properties such as mass and 
spectral functions |l|J^. A large enhancement of low invariant-mass dileptons has indeed been 
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observed in HELIOS-3 [Q] and CERES 0 experiments. Many theoretical studies attributed 
this enhancement either to a dropping p-mass in a medium of high baryon-density [Q , or to 
the collisional broadening of hadronic spectral functions 0 . 

At large invariant masses, dileptons are mostly produced from partonic processes in the 
early stage of heavy-ion collisions. For example, dileptons from the Drell-Yan processes 
dominate the large invariant mass region at the CERN-SPS energy. Thus they can provide a 
benchmark for parton luminosity or the centrality of the collisions in which proposed signals 
of the phase transition to the quark-gluon plasma such as J/'ijj suppression [§ and thermal 
dileptons can be studied via the dilepton channel. However, at collider energies such 
as the BNL RHIC and CERN LHC, the Drell-Yan processes no longer dominate dilepton 
productions at reasonably large invariant masses. Dileptons from decays of charm and 


bottom mesons instead dominate IT^. These large invariant-mass dileptons from heavy 
quark decays may also give us information about energy loss of heavy quarks traveling 


through the dense medium |11-13 


As one can naively expect, both partonic and hadronic processes will contribute to dilep¬ 
tons with invariant masses between the above noted regions, i.e., in the intermediate mass 
region (IMR). An enhancement of intermediate-mass dileptons was observed in S' -|- hP 
collisions at HELIOS-3 [§ and was attributed to secondary meson-meson interactions by 
Gale and Li [|T^. NA38/NA50 experiments [|I^p!6| observed similar enhancement in the 
mass region 1.5 GeV/c^ < M^+^- < 2.5 GeV/c^. The dimuon invariant-mass spectrum in 
p + p collisions is htted by a combination of dileptons from Drell-Yan and decays of charmed 
mesons. Assuming a scaling of binary collisions for hard processes, it was found that dimuon 
spectra m. p + A collisions can be well accounted for but results from S + U and especially 
Pb -|- Pb collisions indicate signihcant enhancement ||T6| , p!7| . It is difficult to argue that the 
excess of IMR dileptons comes from enhanced charm production according to previous es¬ 
timates in which it is concluded that the charm production from secondary parton 

scattering is negligible with a realistic estimate of temperature and gluon fugacity at both 
RHIC and LHC energies. Even if a quark-gluon plasma could be formed at the CERN-SPS 
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energy, the temperature and gluon fugacity are probably too low to produce any significant 
amount of charmed quarks. The charm production in a hadronic gas can also be neglected. 

We focus in this paper on another possible scenario in which IMR dileptons can be 
enhanced in the acceptance of NA50 Pb + Pb experiment. We demonstrate that the mod¬ 
ification of the phase space distribution of D-mesons due to final state rescatterings can 
change the shape of the dilepton distribution in phase space. Since most experiments only 
measure a small part of the whole phase space, the redistribution of the dileptons from de¬ 
cays of charmed mesons can give rise to an apparent enhancement of IMR dileptons within 
the limited acceptance of the experiments. 

In the dense medium formed during heavy-ion collisions, interactions among hadrons 
or partons will lead to partial thermalization of the system. Thus, effects of strong final 
state rescatterings are expected, e.g., in the experimentally measured particle spectra in 
transverse momentum 0. These spectra can be parameterized in an exponential form, 
dN/ mTdrriT oc and are characterized by the inverse slope Te/y. It is found 

that Teff increases with the system size PT| , P2| and is roughly a linear function of the 
particle mass. For example, T^ff for protons increases from 148 MeV in pp collisions to 208 
MeV in central SS collisions, and to 289 MeV in central Pb -|- Pb collisions at measured by 


NA44 The increase is often attributed to a collective transverse flow, which essentially 
comes from final state rescatterings among secondary particles Such effect of final state 
rescatterings which causes the broadening of hadronic spectra, could also happen to 
D-mesoias{D^, , D~ etc). 

To obtain a qualitative feature of this effect, we model final state rescatterings in the 
following way. We assume that D-mesons thermalize with their /oca/environment, where the 
degree of thermalization is parameterized by a temperature T. Due to the lack of knowledge 
about interactions between D-mesons and other hadrons, one cannot accurately estimate 
the value of T. We assume that D-meson spectrum is influenced by final state rescatterings 
in the same way as other hadrons. Thus, T^ff for D-mesons in central Pb + Pb collisions is 
enhanced to about the same value as those for protons and 0-mesons (290 MeV). From this 
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Teff, the parameter T can be constrained to be abont 150 MeV in onr model. Note that T 
is the local temperatnre which differs from the definition of T^ff. Later in the discnssion we 
will show the dependence of onr resnlts on T. 

We calcnlate the open charm prodnction in pp collisions with the latest PYTHIA program 
P5| . We nse the MRS D—' parton distribntion fnnctions with intrinsic parton transverse 

momentnm kt of 1 GeV/c, and with rric = 1.3 GeV/c^. The charm pair cross section in pp 


collisions, cXcl’, is normalized to 9.5 ± 2/ib at Eiab = 200 GeV |T^, which then gives 6.8pb 
at Eiab = 158 GeV. For central Pb + Pb collisions at 158 AGeV, the nnmber of charm pairs 
is then given by 


= Tp,p,(0), 


( 1 ) 


where the overlap fnnction Tpfepb(O) is taken to be 30.4/mb 

We take the Peterson fragmentation fnnction to describe the hadronization of charmed 
qnarks to D-mesons: 

1 


D{z) oc 


( 2 ) 


z[l — 1/^ — e/(l — z)]"^ 

where 2 ; = Pn/Pc, and e = 0.06. D-mesons are then assnmed to be thermalized in its local 
frame with a temperatnre parameter T as the following. For a D-meson with an initial 
fonr-momentnm pD,a = {Ed, Pd) in the Lab frame, hrst a local random momentnm in the 
D-meson rest frame is generated from the thermal distribntion dN/d?p^j^ oc e Then 

the local fonr-momentnm {E\y,p^j^) is boosted back to the Lab frame, i.e., to the frame 


with 7 p) = 1/y 1 — /3p, where (3d = Pd/Ed- This will simulate a locally thermal D-meson 
spectrum with some collective flow in both transverse and longitudinal directions. The local 
temperature T is taken to be 150 MeV if not specihed otherwise in this paper. 

Shown in Fig.|^are the mp spectra within —0.5 < ud < 0.5 and the rapidity distribution 
for initial and hnal D-mesons. All hgures are normalized to one central Pb + Pb event at 
Eiab = 158 AGeV. As expected, initial and final D-mesons have roughly the same rapidity 
spectra, and the inverse slope of the mp spectra at small values of (mp — mp) increases 


4 








significantly after the inclusion of the simulated hnal state rescatterings. The inverse slope 
parameter, T^ff, corresponding to a £t in the region {rriT — G (0,1) GeV/c^ increases 
from 160 MeV to 290 MeV. To show the correlation of charm and anti-charm particles in 
a pair, we plot the azimuthal angle difference, the rapidity difference between the two and 
the correlation in cos Ocs in Fig.g, where cos 6cs is the angle of one D meson relative to 
the beam axis in the rest frame of the charm pair. The correlation in azimuthal angles is 
reduced by random scatterings, however, the rapidity gap between the two particles in a 
charmed meson pair remains almost the same. The distribution in cos 6cs would be flat if 
charm mesons would follow a static-hreball distribution so that there would be no correlation 
between two charm mesons in a pair. However, the charmed meson pair are still strongly 
correlated in rapidity with a big rapidity gap in our model due to the simulated longitudinal 
flow. This is the reason why there is still a strong correlation in cos 6cs despite of hnal state 
rescatterings. In fact, the correlation in cos 6cs still remains strong even if one removes the 
initial correlation in azimuthal angles, 0 i 2 , of the charmed meson pair. 

For semileptonic decays of D-mesons to muons, the average branching ratio of c — 
is taken to be 12%. The muon energy spectrum from H-meson decays is consistent with the 


measurement of the MARK-III experiment 1^^. Since the number of charm pairs, i 


IS 


much less than 1, the dilepton spectrum from decays of uncorrelated charm pairs is negligible. 
In addition, we note that like-sign subtraction is applied in the NA50 data analysis. 

The dimuon invariant-mass, pair rapidity, and which is the angle of one muon 
relative to the beam axis in the dimuon rest frame(Collins-Soper reference frame), are dehned 


as: 


= yipa +p'i y 
Fk+u- = tanh"^ 


cosO^g = 


''z^pair 

pf - A 


z,pair 


Ei^ 


- Plpair\jMl+^-/4: - ml 
where p^^pair = {pf +pf)/{E>^'^ -F E^“). 

Invariant-mass spectra of dimuons from charm pair decays are plotted in Fig.^. The total 
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dimuon yield does not change since the total charm yield does not change. However, there 
are much more high invariant-mass dimuons after hnal state rescatterings. The dimuon pair 
rapidity distribution has little change. 

To consider the effect of the limited acceptance of an experiment, we approximate the 
NASO acceptance by taking the following cuts on dimuon variables: 

> 1.5 GeV/c' e (0.,1.) ,cos0^; G (-0.5, 0.5) , 

and cuts on single muons: 
e^^ €(0.037,0.108), 

E>^> 8. + 1600 X - 0.065)2, ^ (0.037,0.065), 

>8., €(0.065,0.090), 

> 8. + 1300 X (0^ - 0.090)2, ^ (0.090, 0.108), 

where - 2.9 for EL^b = 158 AGeV. 

The invariant-mass, pair rapidity and cos^^^ spectra inside this NASO acceptance are 
shown in Fig. ^ Due to the broadening of the D-meson rriT spectrum, the dimuon yield 
from D-meson decays in the above acceptance is enhanced signihcantly. The enhancement is 
stronger at larger dimuon invariant masses, reflecting the relative change of the rriT spectrum. 
However, the shapes of the pair rapidity spectrum and the spectrum stay roughly 

the same. We dehne the enhancement factor of the dimuon yield from open charm decays 
inside the acceptance, R, as the ratio between accepted dimuon yields after and before hnal 
state rescatterings. It is equivalent to the area under the solid curves divided by the area 
under the dashed curves in Fig. |^b and c. With a value of T = 150 MeV for hnal state 
rescatterings which corresponds to a D-meson rriT inverse slope of 290 MeV, one obtains an 
enhancement factor about 72 = 3.0. 

As we have demonstrated so far, the enhancement factor of IMR dimuons from charmed 
meson decays, R, depends on the D-meson transverse momentum spectrum characterized 
by Te//- This inverse slope is related to the parameter T, which in our model controls the 
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strength of final state rescatterings. As one can easily imagine that the degree of hnal state 
rescatterings must depend on the size of the dense matter in heavy-ion collisions which in 
turn depends on the impact parameter. In principle one can study the effect of rescatterings 
in a cascade model and study the dependence of the effective Tg// or the enhancement factor 
R on the centrality of heavy-ion collisions. To demonstrate such possible dependence, we 
simply vary the effective local temperature T here and study how the effective IMR dilepton 
spectra will change. Show in Fig.|^ are the hnal D-meson rriT spectra for different values of 
temperature T. The inverse slopes, Teff, from these spectra can also be read off in Fig. P, 
where the relations R{T(.ff) and T{Teff) are shown. The two boxes on the left in Fig. prefer 
to the expected dimuon yield by simply scaling up pp results without consideration of hnal 
state rescatterings, thus T = 0 and i? = 1 by dehnition. The two boxes on the right refer to 
the enhanced dimuon yield due to hnal state rescatterings with a default value of T = 150 
MeV. Both the R{Teff) and T{Teff) relation are roughly linear, thus the R{T) relation is 
also almost linear. However, the relation between T or T^ff and the impact parameter could 
be more complex, which we will not address in this paper. 

We would like to point out that the absolute value of the total charm cross section, 
is not important for our study of the enhancement factor, R, for dimuons from charmed 
meson decays; because we normalize i? to 1 for pA collisions and thus acE cancels. We take 


acE to be 6.8/ib just to be consistent with the NA38 data 


Besides hnal state rescatterings, there is also some ehect from initial multiple scatterings, 
which ehectively increases the parton intrinsic kt- By modeling initial multiple scatterings 
in a random-walk picture where the additional average is proportional to the system size, 
we checked that this ehect gives a dimuon enhancement less than 10%. Anti-shadowing 
could be a source for the enhancement too, because at SPS energies the charm production 
probes the region with Bjorken x ~ 0.2, where the gluon anti-shadowing could be close to 
its maximum value |^. However, this ehect is also checked to be no more than 20%. 

In summary, the ehect of hnal state rescatterings on D-mesons in Pb + Pb collisions at 
SPS energies are studied. We found that the dimuon yield from open charm decays inside the 
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NASO acceptance is significantly enhanced as a result. Our model for the effects of hnal state 
rescatterings is very schematic. However, the essential point is the broadening of D-meson 
rriT spectra, which can be roughly represented by the increase of Teff, the inverse slope of 
the D-meson rriT spectrum. More studies can be done beyond this schematic model. The 
centrality dependence of the enhancement can be studied and compared to experimental data 
at various values of total Et- The p± dependence of the enhancement can also be studied by 
applying various p± cuts, and a larger enhancement is expected when more stringent energy 
cut is applied to single muons. However, the ultimate test of this model should come from 
direct measurements of D-meson spectra. 
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FIGURES 



FIG. 1. (a) rriT spectra of initial and final D-mesons at mid-rapidity, (b) Rapidity spectra of 

initial and final D-mesons in the CMS frame. The dashed curve represents initial D-mesons after 
the Peterson fragmentation, and the solid curve represents final D-mesons after the thermalization 
with T = 150 MeV. 
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FIG. 2. Spectra (a) as a function of the difference between the two azimuthal angles; (b) 
as a function of yu, the difference between the two rapidities; and (c) as a function of cosOcSi of 
a charmed meson pair. 
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FIG. 3. Invariant-mass spectra of dimuons from charm pair decays. The dashed and solid 
curves come from D-meson decays without and with final state rescatterings, respectively. 
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FIG. 4. 


(a) Invariant-mass, (b) pair rapidity and (c) cosO^g spectra of dimuons from charm 


pair decays inside the approximate NASO acceptance. 
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FIG. 5. rriT spectra of D-mesons at mid-rapidity after different final state rescatterings. 
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FIG. 6. The dimuon enhancement factor, R, and the parameter T, as a function of the D-meson 
inverse slope, T^ff. The two boxes on the left refer to the case without final state rescatterings, 
while the two boxes on the right refer to the default case of final state rescatterings with T = 150 
MeV. 
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